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ABSTRACT 

A new XMM observation of the outer regions of the galaxy group NGC 5044 indicates hot gas iron 
abundances of only Z-p / Z Q ~ 0.15 between r — 0.2 — 0.4 r v ; r . While the total baryon mass within the 
virial radius may be close to the cosmic mean value observed in rich clusters, the ratio of total iron 
mass to optical light in NGC 5044 is about 3 times lower than that in rich clusters. The remarkably 
low iron abundance over a large volume of the intergroup gas in the outer regions of NGC 5044 cannot 
be easily understood in terms of the outflow of enriched gas in a group wind during its early history 
or by the long term enrichment by the group member galaxies that currently occupy this region. It 
is possible that the stars in NGC 5044 did not produce iron with the same efficiency as in clusters, 
or the iron resides in non-luminous clouds or stars, or the entropy of the iron-enriched gas created in 
early galactic starburst winds was too high to penetrate the group gas of lower entropy. 
Subject headings: X-rays: galaxies: clusters - galaxies: halos - galaxies: formation - cooling flows - 
galaxies: individual: NGC 5044 



1. INTRODUCTION 

One of the most troublesome and least-discussed in- 
consistencies with strict hierarchical cosmological scaling 
is the systematic difference of metallicity between rich 
clusters and smaller galaxy groups from which they pre- 
sumably formed. Most heavier elements are contained in 
the hot intracluster gas and can be measured via X-ray 
spectroscopy. Although these abundance differences may 
require modification of IMF- averaged supernova yields, 
the sense of the group-cluster discrepancy is tha t groups 
conta in less metals per Lb than clusters (e.g., iRenzinil 
1997). Similar conclusions were found from gasdynam i- 
cal models of groups (e.g.. iBrjghenti fc Mathew3 fl9991. 

The ratio of iron mass to optical light in rich clusters of 
galaxies, which are expected to have retained all the met- 
als produced by star s, is Tf« ya , s = My^.^/Ln ~ 0.01 5 
in solar units (e.g., IRenzinil Il997t iLoewensteinl 12004 L 
Adding the iron locked in stars, this increases to Tp e ~ 
0.02. Here we assume an average stellar mass to light 
ratio Tb = (M*/Lb) — 4 and an average stellar iron 
abundance (Zp e ) = 0.7 Zp e ^Q. This value implies that 
about .0035 M^ of iron are pr oduced per M Q of stars 
formed (Kawata & Gibson 2003). The ratio of iron mass 
to optical light (Tpe.gas) for groups is typically le ss than 
about half that of rich clusters fe.g.. lRenzinilll99 7). Even 
with the superior capabilities of Chandra and XMM, X- 
ray observations of bright groups have yet to extend to 
the virial radius, so uncertain extrapolations are required 
to estimate their total iron and baryon content. 

The strong radial iron abundance gradients observed in 
the hot gas of many isolated galaxy groups are particu- 
larly interesting since they suggest that this gas has been 
relatively undisturbed and unmixed by merging. Gas in 
the central parts of groups is enriched primarily by Type 
la supernovae (SNIa) in central E or cD galaxies. How- 
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ever, more distant gas is expected to be enriched by early 
galactic winds from non-central, low luminosity group 
member galaxies. This raises the interesting question 
whether or not the smaller amounts of iron observed in 
hot gas near the periphery of groups can be understood 
as enrichment by SNII and SNIa in non-central galaxies. 

There is general agreement from both ROSAT and 
ASCA that Z-p P , < lZ^ outside the cores o f groups (e.g., 
Finoguenov & Ponman 1999; Buotc 2000a). But because 
of limitations in sensitivity and spectral or spatial resolu- 
tion, these results were considered tentative. Typically, 
ROSAT and ASCA fo und ZyJZp = 1 - 0.5 out to 
r ~ 100 kpc. Recently, iSun et aiT l|2*003') have reported 
iron abundances out to 200 kpc in the galaxy group NGC 
1550. Inspection of t heir figure 4 reveals (resca led to the 
solar abundances of iGrevesse fc Sauva]lll99r1) that be- 
tween 100-200 kpc Z-pe/ Zq m 0.2 but apparently with 
errors that allow for values as high as Z-p a /Z Q = 0.3. 

Previously we reported abundances within the cen- 
tral region of the X -ray luminous group NGC 5044 
ijBuote et al.l l2003blat hereafter Paper 1 and Paper 2 
respectively) finding Z-p e / ' Zq « 1 within r w 30 kpc 
and Z-pe/ Zq ph 0.3 near r w 100 kpc. In this Letter 
we report accurate iron abundance measurements out 
to r k 330 kpc (~ 0.4r V j r ) obtained from a new XMM 
observation of NGC 5044 offset from the center of the 
group. We demonstrate that the very sub-solar values of 
Zp e we obtain at large radius are in serious conflict with 
metal enrichment models for galaxy groups, even in the 
presence of group winds that eject baryons and heavy 
elements from the system. Full details of the observa- 
tions, enrichment models, and derived mass distributions 
of stars, gas, and dark matter will appear in forthcoming 
papers (Papers 3 and 4). We assume a distance of 33 
Mpc to NGC 5044 for which V = 9.6 kpc. 

2. IRON ABUNDANCE FROM XMM 

The 22.5' observation offset due south of NGC 5044 
(Figure P) was performed with the EPIC MOS and PN 
CCD cameras for a nominal 40 ks exposure during AO- 
2 as part of the XMM Guest Observer program. The 
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Fig. 1.— (Left) EPIC MOS images (0.5-5 keV) for the AO-1 (top) and AO-2 (bottom) observations of NGC 5044. Overlaid are the 
annuli used to extract the spectra discussed in this paper. (Right) The 0.5-5 keV spectra of the MOS (lower curve - MOSl:black, MOS2:red 
) and pn (upper curve, green) data in bin #3 fitted with a single APEC thermal plasma. 

TABLE 1 

Temperature and Iron Abundance of the Hot Gas 



| 2D 3D 

n n '"out T Zyc i T Zp c 

Bin (arcmin) (kpc) (arcmin) (kpc) (keV) (solar) (x 2 /dof) (keV) (solar) (x 2 /dof) 



1 10.0 94.8 13.5 128.0 

2 13.5 128.0 19.0 180.1 

3 19.0 180.1 26.0 246.5 

4 26.0 246.5 34.5 327.1 


1.06 ±0.01 0.22 ±0.01 1132.7/1008 

1.07 ±0.02 0.12 ±0.01 554.4/537 
0.98 ±0.03 0.16 ±0.03 727.5/678 
0.81 ±0.03 0.15 ±0.05 742.9/659 


1.06 ±0.01 0.27 ±0.02 1125.5/1008 
1.22 ±0.06 0.13 ±0.03 552.3/537 
0.98 ±0.03 0.16 ±0.03 727.5/678 


Note. — Results of fitting a single APEC plasma model (with Galactic absorption 


) to the XMM spectra. Quoted errors are 



la. "2D" indicates no deprojection was performed. Radial bin 4 is not included in the deprojection analysis (see text). 



background was very quiescent for the entire observation 
allowing for cleaned exposures near the nominal value of 
40 ks for all CCDs. Although no bright, extended sources 
appear in the offset field, ~ 50 faint point sources were 
detected and masked out of the ensuing analysis. 

Since background emission is a sizable fraction of the 
total flux in the offset region, we modeled the back- 
ground rather than using results from blank fields. We 
followed pre vious studies of th e X-ray background with 
XMM (e.g.. lLumb et al]l2002^ and modeled the instru- 
mental background using out-of field-of-view events for 
each CCD. The cosmic X-ray background (CXB) was 
modeled using two thermal plasma components for the 
hot Galactic halo and a power law for unresolved point 
sources. Our CXB constraints are quite similar to t hose 
reported in previous studies (e.g.. lLumb et al1l20 02). 

We extracted spectra from each XMM detector 
(MOSf, MOS2, PN) with sections of circular annuli as 
shown in Figure H and listed in Table Q The center of 
these annuli coincides with that used Papers f and 2. 
(Note that the portions of the annuli that do not lie on 
a detector are excluded from the analysis.) Radial bin 
1 lies immediately outside the region studied in Papers 
f and 2 but is completely contained within the field of 
the AO-f pointing. Consequently, we include the AO-f 
data for radial bin f in the present analysis and model 
the background similarly to that described above for the 
AO-2 observation. Further details of the data reduction 
and background analysis are discussed in Paper 3. 

The source+background models were fitted simultane- 
ously to the spectra of each XMM detector. The source 
models are APEC thermal plasmas with temperature, 



iron abundance (Zp c ), and normalization as free param- 
eters. We use the new standard solar abundances of 
iGrevesse k, Sauval 1)1998(1 . The source parameters and 
the normalizations of the background components were 
fitted separately in each radial bin. Results using the 
best-fitting background model are quoted in this Letter. 

The parameters of the source models obtained from 
fits to the data without deprojection (i.e., "2D" models) 
are listed in Table ^ A single temperature component 
provides a good description of the data. (The fit to the pn 
is not quite as good as to the MOS due to imperfections 
in the background model. But excluding key residuals, 
such as energies between 0.55-0.65 keV and above 3 keV, 
gives an acceptable fit and provides temperature and Zp e 
consistent with the full-band analysis.) Even at these 
large radii the Fe L shell lines near 1 keV are clearly 
visible (Figure^ and allow for precise constraints on Zp c . 
Over most of the AO-2 field Zp e 0.15Z Q which is much 
less than the (emission-weighted) average value Zp e = 
0.44 ± O.O2Z obtained over 48-96 kpc with the AO-1 
data (Paper 2). We note that these low Zp e values are 
robust to errors in the background model; e.g., for radial 
bin 3, we find that the la error on Zp c increases from 
0.03Z Q to 0.05Z Q when accounting for statistical error 
in the background model parameters. Also, accounting 
for weak Fe, O, an d Ne lines due to sol ar wind charge 
exchange emission l|Snowden et al J 12004) would reduce 
slightly our measured Zp e values, particularly in bin 4 

Since Zp D ca n be u nderestimated according to the 
"Fe bias" (e.g.. !Buotell2000rj|) . we investigated whether 
multi-temperature spectral models and spectral depro- 
jection made an important difference in the fitted Zp e val- 
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TABLE 2 
Integrated Masses 



r 


Afgas 


Aagtars 


Mgrav 


A^Fo.gas 


(kpc) 


(lO n M ) 


(1O U M ) 


(1O 13 M ) 


(1O 8 M ) 


327 


12.6 ±0.08 


5.8 


2.00 ±0.09 


3.5 ±0.2 


870 


45.5 


7.4 


3.9 


9.0 



Note. — Errors are 1<t. Quantities listed within r v ; r = 
870 kpc are extrapolations using data within r = 327 kpc. 
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Fig. 2.— Gas iron abundances in NGC 5044. New XMM 
observations for 95 < r < 330 kpc are shown with (open circles) 
and without (triangles) deprojection. The estimated on-the-spot 
total abundance (solid line), is the sum of the gas enrichment due 
to the central E galaxy (dotted line) and the ensemble of non- 
central member galaxies (short dashed line). The upper solid line 
shows Zp c of a typical gasdynamical model with a supernova iron 
enrichment scenario similar to that of rich clusters. The long- 
dashed line shows the on-the-spot abundance when the gas is not 
enriched by SNIa in the non-central galaxies. 

ues. Multi-temperature models (e.g., two-temperature) 
do not improve the fits significantly and do not yield 
larger Zp c values. Deprojecting the data (i.e., "3D" mod- 
els - see Paper 1 for details) under the assumption of 
spherical symmetry also does not significantly improve 
the fits (see Tabled, though the value of Zp c is raised 
by 0.05Zq in radial bin 1. The assumption of spheri- 
cal symmetry allows us to extrapolate the spectra be- 
yond the areas of the CCDs shown in Figure ^ to the 
entire azimuthal range. Such extrapolation magnifies az- 
imuthal fluctuations in the source and background in- 
tensity within the sections represented by the AO-2 field 
and becomes increasingly more important in the outer 
radial bins: the area extrapolation factors are 5 — 6 for 
bins 1-3 and « 9 for bin 4. When bin 4 is included in 
the deprojection, the azimuthally averaged gas density 
does not decrease monotonically with increasing radius. 
Consequently, we exclude bin 4 from the deprojection 
analysis in this Letter. 

These measurements currently represent the most ac- 
curate constraints on Z-p e at radii 200-300 kpc from the 
center of any T ~ 1 keV galaxy group and demonstrate 
Zp c w 0.15Z© at these large radii. Below we use these 
measurements to constrain the global chemical enrich- 
ment of NGC 5044. 

3. INVENTORY OF MASS COMPONENTS 

Estimates of the expected hot gas iron abundance 
require an inventory of the various mass components: 
hot gas, stars, dark matter and iron in the hot gas. 



The integrated masses listed in Table El refer to single- 
phase models of the hot gas except for Mpcgas- Since 
two-temperature models are superior fits to XMM and 
Chandra data within r < 30 kpc and yield larger Zp 
values, the values of Mpcgas in Table El reflect two- 
temperature models in the core of the group and single- 
temperature models elsewhere. However, since the gas 
mass is weighted heavily at large radius, the value of 
-^Fe.gas obtained using a single-temperature throughout 
is quite similar: M Fc = 3.2 ± 0.2 x 10 8 M o . 

We determined the gravitating mass assuming hydro- 
static e quilibrium cl o sely f ollowi ng the procedu res out- 
lined in iLewis et all l|2003ft and iBuote fc Lewis! j.2004). 
The gravitating ma ss is well described by an NFW 
ijNavarro et aljfl997|) model with concentration parame- 
ter, c = 11.1±0.4, and virial radius, r v i r = 870±0.02 kpc. 
(Note: here r vir is defined within an overdensity of 103 
appropriate for the ACDM concordance model.) The 
value for concentration is within reasonable sca tter ex- 
pected for CDM halos (e.g.. iBullock et alJl200lT) . 

To estimate the iron produced by supernovae, we re- 
quire an estimate of the radial profiles of stellar mass 
and density in the NGC 5044 group. The central el- 
lipti cal galaxy can be fit with a de Vaucouleurs pro- 
file ijGoudfrooii et all 11994(1 with effective radius R e — 
62.81" = 10.05 k pc. With RGC3 m agnitudes and extinc- 
tion A B = 0.30 l|Tonrv et al1l200lD we find L B /L B , G = 
4.50 x 10 10 . This corresponds to a total stellar mass 
M^e = T b ,eL b = 3.4 x 10 11 M© where T B ,E = 
Ms, I Lb = 7.5 is ap propriate for a galaxy of this Lb 
l|Gerhard et alJl2001l) using H = 70 km sec -1 Mpc" 1 . 

The central E galaxy is surrounded by ~ 160 smaller 
galaxies (—19 < Mb < —13) that we shall collectively re- 
fer to as "dwarfs" . Only a few dwarf galaxies have known 
radial vel ocities. Of the 163 group me mber galaxies tab- 
ulated bv iFereuson fc S andaee (1990, hereafter FS) we 
have eliminated three (FS 68, 102 & 137) because their 
velocities differ by more than ±800 km s _1 from NGC 
5044. Luminosities Lb have been determined from the 
magnitudes from FS (with Ab — 0.30). According to FS, 
most dwarf members are early types fdE, dSph, SO, etc. ) 
for which we assume (B - V) « 0.7 ijKarick et alJl2003^ . 
Luminosities are converted to mass with T B dF, ~ 4.5, 
typical for early type dwarfs (Gc ha et alJ 12002ft and 
^B.di ~ 2.3 for late types. The total luminosity and 
mass of all non-central galaxies are Lb 4 ~ 1.0 x 10 11 
Lb,q and M* ]£ j w 3.9 x 10 11 Mq. The radial mass dis- 
tribution of the dwarf galaxies is fit with a King profile 
P*A = Po[l + (r/r c ) 2 ]~ 3 / 2 extending out to r t — 570 kpc 
with po = 1.17 x 10 4 M Q kpc" 3 and r c = 130 kpc. Fi- 
nally, using parameterized fits to the integrated masses 
within 327 kpc we extrapolated the masses of each com- 
ponent to r V j r k 870 kpc and list the results in Table El 
For MFe, g as we used a power-law extrapolation based on 
values in radial bins 1-3. 

4. BARYON AND IRON CONTAINMENT 

From Table El the total baryonic mass within the max- 
imum radius observed, r b s w 327 kpc, is Mi,(r & s ) = 
M gas (r obs ) + M*(r obs ) « 1.8 x 10 12 M . The ratio of 
baryonic to total mass at r j }s is therefore fb(r bs) ~ 0.09. 
The baryon fraction /b(r V i r ) « 0.14 extrapolated to the 
virial radius is less certain but is about 85 percent of the 
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cosmic value VL b /Q, m = 0.044/0.27 = 0.16. Only about 
15 percent of the baryons have been lost from the NGC 
5044 group, possibly by group winds at early times if the 
SNII heating efficiency is sufficiently high. 

The total iron mass deficiency of NGC 5044 is much 
larger. The total amount of stellar iron in the central 
galaxy is M Fe .^ E « (Z FeiB )(M* >E /lA) « 2.9 x 10 8 M 
where we assume {Z Fe ^ F ) ~ 0.7Z Fej Q. The total iron 
mass in the dwarfs, M Fe » <j, i s only 1/2 of this, a ssuming 
{Z FeA ) = 0.24 Z Fe& (e.g.. iRakos eFaP 12001^ . With 
these assumptions, the total mass of iron in both stars 
and gas within r f, s and r vlI is M Fe (r i, s ) rj 6.5 x 10 8 Mq 
and M Fe (r viT ) w 13.3 x 10 8 M & . 

If iron was initially produced in NGC 5044 with the 
same efficiency as in rich clusters, we expect that the 
total mass of iron from all past SNII and SNIa within r v i r 
to be M Fe , exp w 0.0035M*(r vir ) « 2.6 x 10 9 M , which 
exceeds the observed iron mass by at least a factor of ~ 2. 
The total iron mass to light ratio, T Fe 0.009, is about 
2 times less than observed in rich clusters. Evidently, 
about 50 percent of the total iron is unaccounted for in 
NGC 5044. The amount of missing iron is much greater 
than the global depletion of baryons, a result that is not 
very sensitive to our extrapolation to r v ; r . 

Figure illustrates the hot gas iron abundance mea- 
surements for NGC 5044 within the 100 - 300 kpc annu- 
lus of the offset region. The iron abundance, Z Fe /Z Q w 
0.1 — 0.2, is significantly lower than the typical abundance 
in rich clusters, Z F£ /Zq ~ 0.4. Is the low Z Fe in this 
distant gas in NGC 5044 consistent with iron-enriched 
outflows from the dwarf galaxies in this same region? 
For an approximate answer to this question we assume 
that the local gas was enriched by dwarf galaxy winds 
according to an "on-the-spot" scenario, where the dwarf 
galaxies (and distant stars of the central E galaxy) are 
surrounded by the same ambient gas that accompanied 
them when they first entered the group halo. Orbital 
models for dwarf member galaxies show that their time- 
averaged radial positions lie near the virial radius at the 
time when they first joined the group. In the on-the- 
spot approximation, the iron abundance in the hot gas 
is found by subtracting the iron in stars from the total 
iron produced by all supernovae over time, 



J Fe.ots 
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where the two terms represent the iron produced by the 
central E galaxy [i = E) and the surrounding cloud of 



dwarf galaxies (i — d). The stellar iron abundance in the 
central E galaxy is app roximated with Z Fe ^ tF (r)/Z Q rj 
0.675 (r fepc /10.05)-°' 21 l|Koha,va,shi fe ArimotollT99l al- 
though the central galaxy does not appreciably enrich 
the offset region 100 < r < 300 kpc. 

The on-the-spot iron abundance Z Fe ^ ots (r) in the offset 
region plotted in Figure 2 is seen to exceed the observed 
abundances by factors of 3 - 4. Indeed this discrep- 
ancy may be underestimated since on-the-spot enrich- 
ment does not include iron expelled from former massive 
group member galaxies (brighter than the most luminous 
surviving non-central galaxy, Mb = —19) that were orig- 
inally in the offset region but have since merged with the 
central E galaxy by dynamical friction. The iron abun- 
dance observed in r > 100 kpc is also very much less than 
that predicted by standard ID gasdynamical models of 
early group w inds in which metals are un iformly mixed 
with the ICM l|Brighenti fe Mathewsf l999). The iron en- 
richment expected in a typical gasdynamical model us- 
ing rich cluster SNII rates shown in Figure 2 results in 
very large ^Fe,gas 

beyond 100 kpc. This is due to (1) 
the outward advection of supernova products from the 
massive central galaxy and (2) the increased efficiency of 
increasing ^Fc.gas by SNIa enrichment when the local gas 
density is lowered by a global group wind. Finally, accre- 
tion of (possibly iron-poor) primordial gas cannot easily 
account for the low observed Z Fe ^ gas . Cosmologically 
accreted gas, which is included in our dynamical mod- 
els, flows differentially inward past the dwarf galaxies by 
less than 70 kpc over 10 10 years. Shocked gas currently 
arriving at r v i r cannot mix deeper into the group IGM 
because of its higher entropy. 

Both the on-the-spot approximation and the detailed 
gasdynamical calculations over-predict the iron abun- 
dance in the offset region. We offer several possible ex- 
planations: (1) the stars in NGC 5044 did not produce 
iron with the same efficiency as in clusters; in particular, 
Type la supernovae in non-central galaxies do not occur 
at the expected rate or failed to enrich the local inter- 
galactic gas, (2) iron has been selectively ejected from the 
group or resides in non-luminous clouds or stars, or, (3) 
high entropy gas enriched and heated by early SNII and 
SNIa explosions may not have penetrated deeply inside 
the NGC 5044 group because of its inherent buoyancy. 
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